Stochastic theory of synchronization transitions in extended systems 
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We propose a general Langevin equation describing the universal properties of synchronization 
transitions in extended systems. By means of theoretical arguments and numerical simulations we 
show that the proposed equation exhibits, depending on parameter values, either: i) a continuous 
transition in the bounded Kardar-Parisi-Zhang universality class, with a zero largest Lyapunov 
exponent at the critical point; ii) a continuous transition in the directed percolation class, with 
a negative Lyapunov exponent, or iii) a discontinuous transition (that is argued to be possibly 
just a transient effect). Cases ii) and iii) exhibit coexistence of synchronized and unsynchronized 
phases in a broad (fuzzy) region. This phenomenology reproduces almost all the reported features of 
synchronization transitions of coupled map lattices and other models, providing a unified theoretical 
framework for the analysis of synchronization transitions in extended systems. 

PACS numbers: 05.45.-a, 05.45.Xt 
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Arrays of coupled oscillators appear ubiquitously in 
Nature, and their mutual synchronization is undoubt- 
edly one of the most intriguing phenomena appearing in 
complex systems Q. Some well-known instances of syn- 
chronizingsystems are chemical reactions 0, neuronal 
networks |3, flashing fireflies 0], Josephson junctions 
electronic circuits |j, and semiconductor lasers Q. In 
particular, coupled map lattices (CMLs) 0, initially in- 
troduced as simple models of spatio-temporal chaos 01 > 
have received a great deal of attention, as models of syn- 
chronization in spatially extended systems (another pos- 
sibility is to analyze discrete cellular automata |iA Ho|h 
In this context, it was observed that distant patches of 
a given CML can oscillate in phase ^jji a phenomenon 
related to the Kardar-Parisi-Zhang (KPZ) equation [J^ 
describing the roughening of growing interfaces . Af- 
terwards, it was also realized that when two different 
replicas of a same CML are locally coupled, they become 
synchronized if the coupling strength is large enough 
B HQ Also globally coupled CML can achieve mu- 
tual synchronization |15| with interesting implications in 
neuro-science 0|. Last but not least, different replicas 
of a CML can be synchronized if they are coupled to a 
sufficiently large common external random noise, even if 
they are not directly coupled to each other [ij, LLa UJ( ■ 

In all these examples, there is a transition from a 
chaotic or unsynchronized phase in which perturbations 
grow, and two replicas evolve independently, to a syn- 
chronized phase in which memory of the initial difference 
is asymptotically lost and replicas synchronize. This syn- 
chronization transition (ST) resembles very much other 
non-equilibrium critical phenomena, particularly transi- 
tions into absorbing states [2(J , and the determination of 
its universal properties has become the subject of many 
recent studies ^J, lltA 1221 l23| as well as the main moti- 
vation of this Letter. Within this context, a major con- 
tribution was made by Pikovsky and Kurths (PK) who 



proposed a stochastic model for the dynamics of pertur- 
bations in locally coupled CML, in which, under very gen- 
eral conditions, the difference-between-two-replicas field 
can be described by the so-called multiplicative noise 
(MN) Langevin equation [22|: 

dt<f>(x, t) = -a<f> - b(j> 2 + DW 2 (j){x, t) + <r<p(x, t)rj, (1) 

where <j>(x, t) is the difference field (or "synchronization 
error"), a, b, a 2 , and D > are parameters, and n a 
delta correlated Gaussian white noise with zero average. 
It is worth remarking that the non-linear saturation term 
was added in by hand, in order to ensure that the 
synchronization error remains bounded from above (i.e. 
of the order of the original field variables) |U HH • 

Eq. Q has been extensively studied, unveiling a rich 
phenomen olog y, including a novel non-equilibrium phase 
transition |24j. This equation is related to: i) the prob- 
lem of directed polymers in random media [lH . \'2'^ and, 
ii) the KPZ equation in the presence of a limiting upper 
wall [2^ |. To verify this last, particularly illuminating re- 
lation, it is useful to perform a Cole-Hopf transformation 
(f>(x,t) — exp[h(x,t)] [23l l24j . obtaining 

d t h{x, t) = -a-b cxph + D\I 2 h + D{Vh) 2 + an, (2) 

which is indeed a KPZ equation plus a bounding term — 
the upper wall — pushing positive values of the height h 
towards negative values [22I |24| . The synchronized or 
absorbing phase corresponds in this language to a de- 
pinned interface, falling continuously towards minus in- 
finity, governed by KPZ dynamics. The unsynchronized 
or active phase, on the other hand, translates into a 
pinned-to-the-wall interface. 

Using the PK theory, it could be naively concluded 
that all ST belong to the MN universality class [2H . 
On the other hand, it has been known for some time 
that many ST in discrete cellular automata (CA) belong 
generically to the directed percolation (DP) universality 
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class |9|. This difference led Grassberger to conjecture 
that the discreteness of the local field variables is the 
factor that determines a ST universality class: MN for 
continuous variables and DP for discrete ones [tj . Follow- 
ing Grassberger, the key difference lies in the fact that 
while for continuous variables one has 'incomplete death' 
(fluctuations can always generate unsynchronized regions 
from synchronized ones [lfj), discrete systems do not al- 
low such a process ('complete death'). This difference 
between MN and DP has been characterized in Ref. [2l| • 

Some counterexamples, however, have recently shown 
that this conjecture does not apply all the way: in some 
cases, CML (continuous variables) can exhibit DP behav- 
ior llil \2i% . or present a discontinuity at the transition 
point [13 • Based on these observations a new global pic- 
ture has emerged: there are two types of ST depending 
on whether the largest Lyapunov exponent, A, is zero or 
negative at the transition. In the first case, the transi- 
tion is MN-like and controlled by linear stability effects 
(A changes sign at the transition). In the second case the 
transition is either DP or discontinuous, depending on 
microscopic details, and the transition is located where 
the propagation velocity of the non-linear perturbations 
changes sign, while A is negative. 

The aim of this Letter is to construct a stochastic 
model, including Eq. P as a particular case, able to 
reproduce in a unified framework all the aforementioned 
phenomenology. With this purpose, we consider the fol- 
lowing Langevin equation 

8t(f>(x,t) = -^p- + DV 2 <t)(x,t)+(t)(x,t)i 1 (x,t), 

0(f) 

V(<t>) = a<t> 2 + b<t) 3 + c<t>\ (3) 

which has appeared in the literature in the context 
of non-equilibrium wetting and depinning transitions 
[25I l2*ij . This is a generalization of Eq. including 
saturation effects in a more general way. In particular, 
for b > 0, the parameter c is irrelevant in the renormal- 
ization group sense, and Eq. behaves as Eq. (|TJ (i.e. 
it belongs to the MN universality class). For b < 0, how- 
ever, new physics emerges |25| (note that in the KPZ 
language b < translates into the presence of an at- 
tractive upper wall). This second possibility is required 
to describe transitions in which a non- vanishing value of 
the synchronization-error field is a local attractor of the 
dynamics. For example, for discrete boolean CA, the 
value 4>(x) — 1 should be a minimum of the local de- 
terministic dynamics, and therefore b < is required to 
describe them within this formalism. Also, for discontin- 
uous cml mm this term describes finite jumps from 
one map-sector to another. 

In Fig ^ we show a plot of the exact single-site effec- 
tive potential in the zero-dimensional limit in the inter- 
face representation (i.e. the stationary solution of the 
associated Fokker-Planck equation after performing the 




FIG. 1: Effective potential in the 0-dimensional case in terms 
of —ft for values of a in the depinned phase. Note the attrac- 
tive effect for b < 0. 



Cole-Hopf transformation; see Ref. [21|). For b < a new 
minimum (absent for b > 0), present even in the synchro- 
nized phase, appears — the KPZ wall becomes attractive. 
As long as & > 0, c is irrelevant [13] j the transition is 
in the MN universality class, and it occurs where the 
depinned interface looses its stability, implying that the 
average velocity (h) changes sign at the transition point. 
It is also well known that (h) coincides exactly with 
the largest Lyapunov exponent, A. Therefore, for b > 0, 
the situation is completely analogous to ST in the MN 
class where A changes sign at the transition. On the 
other hand, for b < |27j the transition can be easily 
shown not to occur at (h) = Instead, there is a 

finite interval of values of a, [a*,a c (b)], in which phase 
coexistence is found, ergodicity is broken, and the final 
state depends upon initial conditions [2^. Within this 
interval, pinned interfaces remain usually pinned, and 
depinned interfaces escape towards minus infinity at a 
finite velocity. Below (above) a = a* (a c (b)) depinned 
(pinned) interfaces lose their stability. Consequently, at 
a c the depinned phase is stable (contrary to what hap- 
pens for b > 0) and (h(a c )) — A is negative, in analogy 
with ST in either the DP class or discontinuous transi- 
tions where A < (as explained afterwards the first-order 
transition could be a transient effect). 

In order for Eq. © to constitute a valid theory for ST, 
a DP critical point should be observed in some param- 
eter range for b < 0. To the best of our knowledge DP 
behavior have not been reported for this equation so far. 
As we will show, DP behavior may emerge for b < as 
the forthcoming numerical results show. 

To simulate Eq. ©, it has been discretized, using the 
Ito prescription, in a one-dimensional lattice, using 1 and 
0.001 as discretization space and time meshes, respec- 
tively. Systems sizes run from L = 50 to L = 2000, 
reaching up to 2 x 10 8 iteration time steps in the longest 
runs. Averages are performed over 100 to 10000 runs 
(up to 2 x 10 6 in some cases). Depending on the values 
explored, different universality classes have been found: 

i) MN: Setting b = c = D = a 2 = 1, we observe a 
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FIG. 2: a) Log- log plot of the order parameter time evolution 
(averaged over surviving, pinned, runs) at the critical point, 
for system sizes (top to bottom) L — 25, 50, 100, 200, 400, and 
1000. Both the main slope, 9 = 0.159(5), and the scaling of 
saturation values, /3/u = 0.245(15) reveal DP-scaling behav- 
ior, b) Log- log plot of the average order parameter as a func- 
tion of a c — a. The full line is a fit to the form (<j>) ~ (a c — a) 13 , 
with an exponent j3 — 0.28(1), also compatible with DP. 



transition around a mll « —0.01 in the MN class |2J|. We 
have also verified that the asymptotic properties do not 
depend on parameter values as long as b > 0. 

ii) First order transition: Changing the sign of b 
(setting it to —3), while keeping the rest of the parame- 
ters as in the previous case, we find a discontinuous tran- 
sition at aft = 1.85(1). In this re gim e we recover all the 
phenomenology reported in Ref. 25]: an apparent first 
order phase transition, a broad region of phase coexis- 
tence, triangular patterns in the interface representation, 
etc. However, by enlarging the value of c the attractive 
well becomes less deep (see Fig. QJ, the discontinuity of 
the transition is reduced, and it becomes more and more 
difficult to clearly establish the order of the transition. 
Indeed, following a claim by Hinrichsen |28j it could well 
be the case that this discontinuous transition is only a 
transient effect, asymptotically crossing-over to DP for 
very large system sizes. 

iii) DP: In order to weaken the mechanism leading to 
a first order transition the mechanism responsible for 
destabilizing depinned patches — we consider a shal- 
low well, c = 4, reduce the diffusion constant to D = 0.1, 
and keep b = —3 and a = 1. In Fig[2K) we show how the 
order parameter decays at criticality (a c = 0.1865(3)) 
with a DP exponent 9 w 0.16 for L = 1000, and also 
how it saturates for smaller system sizes. From the scal- 
ing of the saturation values at criticality, we estimate 
[3/v = 0.245(15), in good agreement with the DP value 
(3/v = 0.252... |2fJ. As an independent check of uni- 
versality, we have also studied the dependence of ((f)) 
on a at a fixed system size L = 1000. A fit to the 
form {(f>) ~ (a c — a)P provides the independent estimates 
a c = 0.1871(5) and (3 = 0.28(1), again compatible with 
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FIG. 3: Time evolution of the moment ratios reported in the 
text for L = 25, 50, 100 and 200 at the critical point. Note the 
convergence for large sizes to the expected DP values (straight 
lines). 

DP, see Fig. \^p). To further confirm DP scaling, we 
have computed different moment and cumulant ratios at 
criticality. It has been recently shown that moment and 
cumulant ratios are universal at the DP fixed point, and 
they are not strongly affected by finite size effects. In 
particular, we have measured r\ — (m.2 — m\)/m\ (anal- 
ogous to Binder's cumulant for these transitions) and 
T2 = m^/rri^, where rrii is the i-th moment. In Fig. |3| 
we show their time evolution at criticality for different 
system sizes. Observe the convergence to the extrap- 
olated DP values for infinite lattices n « 0.173 and 
r 2 ~ 1.554 |2^]. Slight deviations from the values re- 
ported in Ref. |2^| are due to lack of precision in the 
determination of the critical point. 

We have also measured a different order parameter: 
the density of lattice sites with h smaller than a certain 
reference threshold — i. e. the number of sites pinned at 
the attractive wall. We have checked that it scales as (</>) 
(DP) at criticality, independently of the arbitrary thresh- 
old. Also, as a complementary test to verify the genera- 
tion of DP scaling for b < 0, we have studied the nature of 
fluctuations. In particular, we determine the distribution 
of the output order parameter, P{(4>{t + 1))), after one 
iteration as a function of the input order parameter value 
(4>(t)). It is not difficult to verify that its variance is pro- 
portional to y/4>(t), as expected for the DP class [i^ • In 
other words, the presence of an attractive wall generates a 
DP-like noise \/(f>(x, t)i~i, to be included in Eq. While 
in the absence of such a wall the leading fluctuations are 
those intrinsic to the interface (linear in (/>) the attractive 
wall generates a new noise term, related to the fluctua- 
tion of the number of sites pinned by the potential well. 
This noise dominates the asymptotic behavior, driving it 
from MN-type to DP. Let us remark that even though 
the transition is DP-like in this regime, it has some first- 
order flavor. Indeed, as stressed before, there is a finite 
phase coexistence interval [a* , a c ] , within which pinned 
and depinned interfaces coexist with a phenomenology 
very similar to spatio-temporal-intermittency |25j . For 
a < a* (a > a c ) only the pinned (depinned) phase is 



4 



stable. A qualitatively similar ST have been found for 
CMLs and the term fuzzy transition has been coined to 
describe it |l7|. 

From the theoretical side, we have performed a renor- 
malization group perturbative calculation of Eq. © along 
the same lines followed to study Eq.Q |24j. We find a 
weak coupling trivial fixed point (above d = 2), and run- 
away trajectories above some separatrix analogously to 
what happens for Eq.QJ [24| (in particular, in d = 1 
these are the only solutions). The possible existence of a 
DP strong-coupling fixed point cannot therefore be set- 
tled within this perturbative approach. This calls for new 
theoretical studies (see |3fij). 

To sum up, all the observed universal features of 
synchronizing CML are well reproduced by a simple 
Langevin equation, Eq.J2J), including a potential with up 
to quartic terms, diffusion, and multiplicative noise, that 
is equivalent to a KPZ equation with a bounding upper 
wall. A MN transition is found for b > (non-attractive 
wall), while for b < (attractive wall) first order or DP 
transitions are observed depending upon parameter val- 
ues. In the MN case the scaling is controlled by inter- 
face fluctuations, while in the second case it is domi- 
nated by fluctuations of the number of sites pinned by 
the attracting wall (generating DP-type of noise). This 
second mechanism has been shown to exhibit DP-scaling 
(probably, DP with many absorbing states [13)01) or, 
alternatively, to induce (apparently) discontinuous tran- 
sitions; in both cases, the transition is "fuzzy": there is 
a finite coexistence region between synchronized and un- 
synchronized phases. In this way we have constructed 
a general framework to study ST in spatially-extended 
systems and, in particular, in coupled map lattices, suc- 
cessfully reproducing all their universal critical features. 
It would be very interesting to derive analytically the 
Reggeon field theory Langevin equation describing DP 
|20j from Eq. Q in the adequate parameter range, and 
to develop clean-cut criteria establishing the nature of 
the transition as a function of bare parameter values (i. 
e. the full phase diagram) . These issues will be addressed 
in a future publication. 
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